Scanning the critical fluctuations 
— application to the phenomenology of the two-dimensional XY-model — 



Ricardo Paredes VJ'+ and Robert Botet* 
t Instituto Venezolano de Investigaciones Cientifica, Centro de Fisica, 
Laboratorio de Fisica Estadistica Apdo21827 1020 A Caracas, Venezuela 
* Nano Structured Materials, Delft University of Technology, Julianalaan 136, 2628 BL Delft, Netherlands 
* Laboratoire de Physique des Solides Bat. 5 10, CNRS UMR8502 - 
Universite Paris-Sud, Centre d'Orsay, F-91405 Orsay, France 
(Dated: February 6, 2008) 

We show how applying field conjugated to the order parameter, may act as a very precise probe 
to explore the probability distribution function of the order parameter. Using this 'magnetic-field 
scanning' on large-scale numerical simulations of the critical 2D XY-model, we are able to discard 
the conjectured double-exponential form of the large-magnetization asymptote. 
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Introduction. — Derivation of the complete equa- 
tion of state of a many-body system is generally a 
formidable task. When the system may appear under 
various phases at the thermodynamic equilibrium, this 
problem requires knowledge of the exact probability dis- 
tribution function (PDF) of its order parameter. Despite 
a number of attempts, just a few instances are available 
0. Even the exact PDF for the 2D Ising model is still 
unknown. 

Within this context, the critical point is very particu- 
lar, since the universality concept tells us that only a lim- 
ited information is needed to obtain the complete leading 
critical behavior. For instance, general arguments give 
precisely the tail of the critical PDF, P(m), for the large 
values of the order parameter, m, namely [2j: 



form: 



P{m) 



(1) 



with c a positive constant and 6 the magnetic field critical 
exponent, or the distribution of the zeros of the Ising 
partition function in the complex magnetic field 3] (such 
a partition function is Fourier transform of the PDF). 

In the present work, we explain how the real magnetic 
field can be generally used as a very accurate probe to 
scan quantitatively the zero-field PDF tail, exemplifying 
the method with the critical 2D XY-model. By the way, 
we will see that the popular double-exponential approxi- 
mation of the PDF for this system cannot be correct, and 
we provide alternative approximation which is consistent 
with the critical behavior. Consequently, our results dis- 
card possible fundamental connexion between this mag- 
netic model and the field of extremal- values statistics. 



Former approximation of the magnetization 
PDF for the critical 2D XY-model. In a series 
of recent papers 0, 0, IE > it was argued that the PDF 
P(m) of the magnetization m of the 2D XY-model at the 
Berezinskii-Kosterlitz-Thouless (BKT) critical tempera- 
ture, could be approximated by the generalized Gumbel 



P(m) oc exp (b a z a - \ cr e a ' jZa ) 



(2) 



where the reduced magnetization: z a = (m — (m))/o~ is 
used. From low-temperature spin-wave theory and direct 
numerical simulations, one obtains [f|: 



1.105 



1.74 



0.69 



(3) 



It was regularly noticed [5j that the form J5j cannot be 
the exact solution of the corresponding statistical prob- 
lem, even if a number of analytical arguments as well as 
numerical simulations show convincingly that this trial 
function is indeed close to the exact solution. Moreover, 
Eq. © is appealing, as it suggests connexion between the 
critical 2D XY-model and the statistics of extreme vari- 
ables |8|. Therefore, the question of a possible bridge be- 
tween these two active fields of statistical physics should 
be examined precisely. On the other hand, Eq.(|2Jl is in- 
consistent with the general behavior QJ, since 5 = 15 
for the 2D XY-model. The question to know whether 
relation Q is true or wrong for this system, is then fun- 
damentally important. We will examine hereafter these 
two questions. 

Two alternative hypothesis. — We consider the 
2D XY-model |9j on a square lattice of size LxL with pe- 
riodic boundary conditions. The N = I? classical spins 
are confined in the x-y lattice plane, and they interact 
according to the Hamiltionian: H = —JJ2<ij> S< ' Sj, 
where J > is the ferromagnetic coupling constant 
and the sum runs over all nearest-neighbor pairs of 
spins. Eventual critical features are characterized by the 
singular behavior of the scalar magnetization per site: 

m = jf \J Si) 2 , which is a non-negative real number. 
We define also the instantaneous magnetization direc- 
tion as the angle ip such that: £LSf = mN cos yb and 
V\ S? = m N sin ijj. 

There is a continuous line of critical points for any 
temperature below the critical BKT temperature Tbkt 



<m> c P c (m) ln(<m> c P c (m)) 




FIG. 1: PDF of the magnetization for the 2D XY-model 
at the critical temperature Tbkt, plotted in the first-scaling 
form The scaling law is confirmed for L = 64 (stars) and 
L = 128 (circles), while the L = 16 (continuous line) shows 
finite-size deviation. Wolff's single-cluster algorithm was used 
fl4| . Each data set corresponds to average over 25,000,000 
independent realizations. 

[Icj . In this region, < T < Tbkt, the system is critical, 
and asymptotic (i.e. L — > oo) self-similarity results in the 
so-called first-scaling law |l2j: 

(m)P(m) = $t(zi) , with zi = — — , (4) 

(m) 

and $t is a scaling function which depends only on the 
actual temperature T . Under this form, the hyperscaling 
relation, (m)/a = est, is automatically realized. Eq.Q 
is a sequel of the standard finite-size scaling theory 
[l3| . but it is highly advantageous that does not 
require knowledge of any critical exponent. FIG.Q 
gives numerical exemplification of the first-scaling 
law at Tbkt, and illustrates the overall shape of the 
distribution $ c (zi) (hereafter, the index 'c' refers to the 
BKT critical point, T = T B kt)- 

We separate the free energy T of the 2D XY-system 
at equilibrium (temperature T = 1/(3) into the sum of 
a regular part describing the small values of the magne- 
tization, a sing ular part vanishing as the essential 
singularity 16, 17] when T — > Tbkt, and a regular part 
for the large values of the magnetization, namely: 

l3T{m) = ip Q (m/ (m)) + (p s {rn/ (m)) + ipoo^m/ (m)) . (5) 

Clearly, discussion on the system behavior can be carried 
out either through the free energy © or the first-scaling 
law since: In P(m) — — /3.F(TO)+constant term. 

The regular small-m tail. - As the singular behavior 
should vanish at the BKT transition, we study first the 



FIG. 2: Part m < (m) c of the logarithm of the scaled PDF 
@ vs z\, for the 2D XY-model at T = Tbkt- The solid 
straight line is the best fit: ln((m) c P c (m)) = &izf+cst for 
the L = 128, 2i < 0.8, data. Numerically, bi = 12.7. Same 
symbols as in Fig0 

regular small-m behavior of P(m) at T — Tbkt- Numer- 
ical results for P c (m) are shown on FIG[21in the form 
They suggest the leading form: 

lnP c (m) w bi (m/(m) c ) 2 . (6) 

Close to the most probable value of the magnetization, 
Eq.|J2J and Eq.@ are indeed consistent each other as the 
latter writes : lnP c (m) = cst + 2biz a /((m) c /a c ) + 0(z%), 
in which we recognize the term linear in z a . 

The singular small-m tail. We consider now the 
singular part of the free energy through the combination: 

In ((m)P(m)) - In ((m) c P c (m)) (7) 

vs the reduced magnetization z\ = m/(m). The data 
plotted in FIG|21 suggest a cubic ^-behavior: 

« C ( T )4 . (8) 

for every T < Tbkt, and for the values of m smaller than 
the mean. Moreover, c(Tbkt) = 0. 

The large-m tail at the BKT point. - Instead of us- 
ing multicanonical Monte-Carlo simulations 01 which 
need too large system sizes to conclude 0], we consider 
static in-plane magnetic field, H, as a probe to study the 
features of the PDF for the large values of the magneti- 
zation. Indeed, as the intensity of H increases, the most 
probable magnetization, m^-, as well as its mean value, 
(ma), explores larger values of the PDF tail. We con- 
sider two alternative forms for the critical tail, namely: 

• the 'Gumbel-like' shape J2J) - noted below: ' hypoth- 
esis (G)' -, which writes in the first-scaling form: 

* c (*si) ~ exp(-A e Q ° 21 ) for Zl -> oo , (G) 



3 



ln(<m>P(m) /<m> c P c (m)) 




FIG. 3: Part m < (m) of the logarithm of the scaled PDF, 
corrected by the regular part of the free energy (see Jjj), for 
L — 16 and four different temperatures: T = 0.3 (circles), 
T = 0.6 (squares), T = 0.8 (diamonds), and T = 0.885 (stars) 
which is close to the critical temperature (Tbkt ~ 0.893 [2(J). 



The plot is versus z\ 
best fits Eq.©. 



(m/(m)) . The straight lines are the 



with a = a a x ((m) c /a c ) (w 16.4 from © and Table Q}, 
and Ao = \ a e~ a °. It is the form suggested in 

• the 'Weibull-like' critical shape - noted below: ' hy- 
pothesis (Wy -, which is [2l]: 

$ c (zi) -exp^Aiz^ 1 ) for Z\ — > oo , (W) 

with Ai a positive parameter, and 5 + 1 = 16 |2^ . 

Let <f) be the direction of H with respect to the x- 
axis (i.e. H = (H cos 0, sin </>)). According to gen- 
eral thermodynamics, the magnetization PDF is given 
by: P c (m, H) oc exp(— (3 C T + f3 c L 2 mH cos(tp — 0)), with 
the field-less free energy T . Therefore, the most prob- 
able magnetization, m* H , is the solution of the equation 
dP c (m, H) I dm — for a given value of H. As the instan- 
taneous magnetization direction ip should coincide with 
the magnetic field direction <j> for the large systems, we 
use cos(ip — (f>) r* 1. Rewritten in terms of the auxiliary 
variables X = H/(m) s c and Y = H/m* H 5 , Eqs.©,®, 
with hypothesis (G) or (W), result respectively in: 



x n , fx 



l/S 



Aoa oe *°W y ) 1/5 (9) 



= M(S+l)y , (10) 

which arc implicit equations for the most probable mag- 
netization, m* Hl (written in the combination Y) vs the 
magnetic field H and the system size N (written in 
the combination X). The constant A is such that: 



For the large magnetic field, m* H is expected to be 
much larger than (m) c , that is: X/Y ^> 1. Consequently, 
the solution of Eq.|0 is: 



Y = a s X/(\nX + C) s 



(11) 



where C — — ln(AAoao) ~ 14.2 is a positive constant. 
Within the hypothesis (W), one has: (X/Y) 1/S < X/Y, 




FIG. 4: Double-logarithmic plot of H/(iriH) 5 vs the reduced 
magnetic field HL VH with: yu =28/(5+1). These two vari- 
ables are convenient for the numerical simulations and simply 
related to the variables X and Y of the text: H/{mH) S = 
Y x (m* H /(m H )) S and HL VH = X x (Af3 c )- 6/{s+1) , with the 
constant values: (m* H /(m H )) S ~ 1 and (Af3 c )~ S/(S+1) « 0.96. 
The dashed curve is the solution of the first Eq.@ (cor- 
responding to the hypothesis (G)), while the dotted line is 
Eg. 1121 . with Ai = A CT in agreement with 1141 . The system 
size goes from L = 16 up to L = 512. Each point corresponds 
to an average over 100,000 independent realizations |23|. 

such that l(TU|) shows that Y is asymptotically a constant: 



Y^AXi(6 + l) 



(12) 



A~ 



\S+1 



1.07. 



So, increase of Y with the intense magnetic field should 
be interpretated as failure of (W). 

Inference from the numerical simulations. 

Both solutions, (fTT|) and (|P2|l . are drawn on FIGQ] 
in comparison with the results of large-scale numerical 
simulations of the 2D XY-model with the in-plane mag- 
netic field at the BKT temperature. It is clear that the 
numerical simulations are consistent with the hypothesis 
(W), while the double-exponential tail (G) should be dis- 
carded. This suggests the following form of the critical 



4 



TABLE I: Temperature, system size, average magnetization 
per spin, ratio of average magnetization to standard devia- 
tion. The best fit for the latter is: (m) c /a c = 14.81 — 21.5/L, 
at the BKT temperature, T BKT = 0.893. 



T 


L 


(m) 


(m)/a 


0.3 


16 


0.923218 


66.958 


0.6 


16 


0.836307 


29.249 


0.8 


16 


0.764091 


18.260 


0.885 


16 


0.723259 


13.907 


0.893 


16 


0.718814 


13.467 


0.893 


32 


0.662819 


14.119 


0.893 


64 


0.611181 


14.486 


0.893 


96 


0.582217 


14.583 


0.893 


128 


0.563209 


14.644 


0.893 


256 


0.518921 


14.687 


0.893 


512 


0.478045 


14.829 



We chose the critical 2D XY-model as a debated exam- 
ple to treat with this method. Indeed, a former double- 
exponential approximation of the magnetization PDF in 
the 0-magnetic field is found to be inconsistent with the 
critical behavior of the system - though correct near 
the most probable magnetization -. Moreover, this ap- 
proximation being taken from another field of statistical 
physics, could mislead, as it suggests hidden link between 
these two fields. The new proposed approximation cor- 
rects these flaws. 
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PDF for the 2D XY-model: 



P c {m) oc e " l2: i~ Al2 i , zi=m/{m) 



(13) 



Below the BKT critical temperature, additional term 
+c(T)zf should appear in the exponential. 

In order to understand the origin of the approximation 
P)l. let us change the reduced magnetization according 
to: zi = 1 + z a / ((m) / a). At Tbkt-, and for the small 
values of z a / {(m) c / a c ) (recall that (z a ) = 0, and that 
(m) c /(T c w 14.8 is a rather large number), we obtain: 

P c (m) CX e 2b ^-/((™)c/^c)-A 1 (l+ 2 , T /((m) c / ( T c )) 16 _ 

Writing then 1 + z a /((m) c /a c ) w e ^/((m) I! / CTc ) ) yields 
Eq.(0l, provided the following relations are verified: 



1G 



2&i 



{m) c /a c 



{m) c /a c 



\* = Ai , (14) 



So, Eq.J2J appears to be a good approximation around 
the most probable magnetization, but is inconsistent with 
the general critical relation (m# ) oc H 1 ' * unlike Eq. (|13fl • 
By the way, the conjectured relation pj b a ja a = tt/2 
writes simply: b\ = 47r, that we accept here as a new 
conjecture (numerically: b 1 » 12.7, see FIGJ2J- 

Conclusion. — In this Letter, we explained how us- 
ing the field conjugated to the order parameter provides 
unique information about the tail of the probability dis- 
tribution function of the order parameter. This is of ma- 
jor importance for the critical systems, since the shape of 
the tail is directly linked to the value of a critical expo- 
nent. Therefore, this general method provides alternative 
way to calculate or measure the critical exponent S. 
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